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Abstract. For a class of nonconservative hyperbolic systems of partial dif- 
ferential equations endowed with a strictly convex mathematical entropy, we 
formulate the initial value problem by supplementing the equations with a ki- 
netic relation prescribing the rate of entropy dissipation across shock waves. 
Our condition can be regarded as a generalization to nonconservative systems 
of a similar concept introduced by Abeyaratne, Knowles, and Truskinovsky for 
subsonic phase transitions and by LeFloch for nonclassical undercompressive 
shocks to nonlinear hyperbolic systems. The proposed kinetic relation for non- 
conservative systems turns out to be equivalent, for the class of systems under 
consideration at least, to Dal Maso, LeFloch, and Murat's definition based on 
a prescribed family of Lipschitz continuous paths. In agreement with previous 
theories, the kinetic relation should be derived from a phase plane analysis of 
traveling solutions associated with an augmented version of the nonconserva- 
tive system. We illustrate with several examples that nonconservative systems 
arising in the applications fit in our framework, and for a typical model of 
turbulent fluid dynamics, we provide a detailed analysis of the existence and 
properties of traveling waves which yields the corresponding kinetic function. 



1. Introduction 

Certain nonlinear hyperbolic models arising in continuum physics and, especially, 
models describing complex fluid flows do not take the standard form of conservation 
laws but, instead, are nonlinear hyperbolic systems in nonconservative form 

(1.1) dtu + A{u) d^u = Q, xeR, t>0. 

Here, u = u{x, 6 f2 is an unknown field taking values in a (convex and open) 
domain D, C R^, while the matrix- valued field A = A(u) is given and, for each 
state M, admits A'^ real and distinct eigenvalues. It is well-known that nonlinear hy- 
perbolic equations do not admit smooth solutions since propagating discontinuities 
arise in finite time even from smooth initial data. For conservative systems, weak 
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solutions in the sense of distributions are sought for. However, for nonconservative 
systems the distributional definition does not apply. A suitable notion of 

weak solution was proposed Dal Maso, LeFloch, and Murat and the nonlin- 
ear stability of such solutions was investigated therein. Nonconservative hyperbolic 
systems have been the subject of active research in the past fifteen years. The 
theory covers the definition of weak solutions [56l [36l |37l [2TJ |42] , the existence of 
solutions to the Riemann problem [36j [21] the initial value problem [36» iM. i2Qj , the 
uniqueness of bounded variation solutions [SI [JT] , and their approximation via finite 
difference schemes [30[ I15[ 145] . In addition, many nonconservative models arising in 
continuum mechanics have been systematically investigated, as such model play an 
important role in the modeling of multi-phase flows and turbulent fluid dynamics 
[3l[a[8l[9l[l6l[l7]. 

Building upon the above works, our purpose in the present paper is to consider 
a restricted class of nonconservative systems of the form (|l.ip , characterized by the 
property that a large family of additional entropy functions (conservation laws) is 
also available. In other words, the systems to be considered below formally have 
a conservative form if nonlinear combinations of the given equations are allowed. 
However, the physical modeling dictate that nonconservative equations should be 
kept and it is precisely under these conditions that a "kinetic relation" , as we 
propose in the present paper, should enter into play. 

The kinetic relation were initially introduced by LeFloch [40] for hyperbolic sys- 
tems of conservation laws in order handle nonclassical undercompressive shocks, 
following earlier works by Abeyaratne, Knowles [1] , and Truskinovsky [55j for sub- 
sonic phase transitions. See |39l [40l [HI [43] for details. 

The concept of a kinetic relation for nonconservative systems discussed now was 
actually introduced first by the authors in an unpublished manuscript |10| . Later 
on, this concept was investigated numerically by Aubert, Berthon, Chalons, and 
Coquel |3l [9j [T7] , and the control of the numerical dissipation of finite difference 
schemes was extensively addressed. Our purpose in the present paper is to pro- 
vide the required theoretical framework and demonstrate that the kinetic relation 
provides an efficient tool to handle complex fiuids. 

Recall that the design and the properties of difference schemes suitable for the 
numerical approximation of nonconservative systems (jl.ip is very challenging. The 
main source of difficulty lies in the fact that shock waves to nonconservative sys- 
tems are small-scale dependent and the dissipation terms induced by the numerical 
discretization tend to drive the propagation of the shocks. This phenomena was 
rigorously analyzed for scalar equations by Hou and LeFloch 30, . On the other 
hand we emphasize that Glimm scheme and front-tracking algorithms do not con- 
tain any numerical dissipation and, actually, have been proven to converge to the 
correct solutions [36, 44, 41 . The method based on the kinetic relation proposed in 
the present paper allows one to extend to nonconservative systems the conclusions 
made for nonclassical shocks in [27l[45] (and the references cited therein). 

We begin with a general discussion of nonconservative hyperbolic systems arising 
in continuum physics in order to motivate our general approach proposed in the 
forthcoming section and developed on selected examples in the rest of this paper. 
The models of interest here naturally stand in a nonconservative form, and this is 
a direct consequence of simplifying assumptions which are made in the derivation 
of these models; these assumptions are also necessary if a tractable model is to be 
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found. Such assumptions typically originate in averaging procedures that intend to 
bypass the description of intricate mechanisms taking place at microscopic scales. 
The small scale fluctuations that are thought to be of lesser interest induce dis- 
sipative and/or relaxation phenomena at the macroscopic level, and can also be 
accounted for as source-terms. 

Most (if not all) nonconservative hyperbolic models arising in the applications 
admit (several distinct) entropy balance laws which are consistent with the underly- 
ing dissipative and relaxation mechanisms. These additional balance equations, as 
we will show, provide a natural approach to formulating additional generalized jump 
conditions built from entropy rate productions. Moreover, these entropy functions 
are sufficient in number to allow for a complete set of jump relations. 

The objectives and results in this paper are as follows: 

First of all, as mentioned above, we restrict attention to a class of nonconser- 
vative systems (defined in Section |3] below) which encompasses, however, most of 
the models encountered in the applications. To motivate the proposed class of sys- 
tems, we observe that, in the applications we have in mind (e.g. multi-phase and 
multi-fluid models): (1) all but one of the equations can be rewritten in a con- 
servative form and, moreover, (2) the system is endowed with a mathematical 
entropy, i.e. a (strictly convex) nonlinear function U = U{u) corresponding to an 
additional conservation law satisfied by all smooth solutions. 

For such systems, the concept of weak solutions introduced by Dal Maso, LeFloch, 
and Murat [21 can be simplified. Therein, a family of Lipschitz continuous paths 
was necessary to uniquely define the nonconservative product A{u) dxU associated 
with the vector-valued field u. In contrast, for our particular class of nonconser- 
vative systems, one nonconservative product between scalar-valued functions, only, 
must be defined. This structure allows us to simply supplement the model (II. 1[) 
with an additional algebraic scalar equation which, for each shock wave, determines 
the entropy dissipation rate associated with the entropy U. We call this additional 
jump condition a kinetic relation and the entropy dissipation function a kinetic 
function. In Section [3] below, a precise definition is given. The main result of this 
section is a proof of the existence of a solution to the Riemann problem for (jl.ip 
which satisfies the prescribed kinetic relation. Our proof is a generalization of an 
argument given in Dal Maso, LeFloch and Murat [21] in the setting of general 
families of paths. 

It is remarkable that many models of interest arising in the applications take 
the form considered in Section |3] below, and this will be illustrated in the following 
Section [2] In Section |4l we focus on a model of particular importance, which arises 
in turbulent fluid dynamics. Taking into account the dissipation terms induced by 
the physical modeling, the existence and properties of associated traveling waves 
are established. In Section we characterize the right-hand states of traveling waves 
with fixed left-hand states, which leads us to the desired kinetic relation. In turn, 
this provides us with the kinetic function needed to apply the general theory in 
Section |3l 
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2. NONCONSERVATIVE SYSTEMS IN FLUID DYNAMICS 

To show the structure of the nonconservative systems of interest, it is worth to 
begin with the shallow water equations with topography 

dtp + dxipv) = 0, 

(2-1) p2 

dt{pv) + dx{pv^ + g —) ~ gpdxa = 0, 

where p and v are the mass density and the velocity of the fluid, respectively, and 
the (prescribed) topography function a : R — > R+ depends on the spatial variable 
X and is assumed to be solely piecewise Lipschitz continuous. Here, g is the gravity 
constant. The product gp dxO. is nothing but a nonconservative product which is 
not defined in a classical sense at points of discontinuity. 

By setting u :— {p,pv), weak solutions should obey the following entropy in- 
equality 

dtU{u,a) + dxJ'{u,a) < 0, 
(2.2) U{u,a) -.^ pE{v) + pa, e'(p) = 



^3 



T{u, a) :— p — + pe{p) v + p{p)v + pv a. 



(2-3) _ , , _ 2 



Another model with a closely related structure is 

dt{ap) + dx{apv) = 0, 

dt{apv) + dx{apv^ +ap{p)) - p{p) dxa = 0, 

which describes one- dimensional nozzle flows as well as compressible Bows in porous 
media. Again, the function a : R — R is solely piecewise Lipschitz continuous, and 
denotes here the nozzle cross-section or the porosity function, respectively. 

By setting u :— {ap,apv), weak solutions to (|2.3p should obey the entropy in- 
equality 

dtU{u, a) + dxJ'ivi, a) < 0, 

(2-4) U{u,a)^a^p—+ape{p), 

J"(u, a) = {U (u, a) + p{p)) V. 

The systems (12. 1[) and (|2.3|) and closely related models with source-terms have 
received considerable attention over the past decade, from, both, analytical and 
numerical standpoints. We refer the reader to [37] (connection with the theory of 
nonconservative systems), [23l [24l |4l [151 fT2^> 'S?, "3? (approximation by finite dif- 
ference or finite volume schemes), and ^8, 49, 31, 25, 46, 47 (construction of a 
Riemann solver). In particular, we refer the reader to Bouchut [12 for a compre- 
hensive review and to the references therein. 

We observe here that both models (|2.1I) and (12. 3p fall within the following class 
of nonconservative hyperbolic models with singular source-term 

(2.5) dtu + dxi{u,a) - g{u,a)dxa ^ 0, 

where a is a given (piecewise Lipschitz continuous) function of the spatial variable 
X and the unknown map u takes values in a convex and open domain il^ C R^, 
while f : fiu x R — > R''^ and g : flu x R R^ are given smooth mappings. 
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Motivated by the structure of the above two examples, especially the entropy 
inequalities (|2.2p and (|2.4|) and in order to develop a general theory we assume 
that the hyperbolic system (I2.5P is endowed with a (sufBciently smooth) entropy 
function U : 0,^ x M ^ M and a corresponding entropy flux J- : il^ x R — )■ R, so 
that solutions to (|2.5p satisfy the entropy inequality 

(2.6) dtU{M,a)+d:,T{\i,a) <Q. 

The principal examples of interest arising in the form ()2.5p in the applications do 
admit such an entropy. 

The above class is known to include, after the seminal work by Greenberg and 
Leroux [23) . the class of hyperbolic systems with source terms: 

(2.7) 9tu + a,f(u) =g(u), 

which, by introducing the (rather trivial function) a{x) = x, indeed take the form 
(cf. (1231) 'I 

(2.8) + (u) - g(u) = 0. 

This nonconservative reformulation is useful for designing "well-balanced schemes" , 
which properly account for the competition between the source term and the differ- 
ential hyperbolic operator in the large time asymptotic t — > +c>o. (See [23 1I241IT2] .) 
The (somewhat artificial but useful) system p.8|) admits a conserved entropy in 
the scalar case n — 1, provided the source g does not vanish, namely it suffices to 
define W(u) := l/g(w) and J^'(u) := f'{u)/g{u). 

As advocated by LeFloch [37j for the equations for nozzle flows (j2.3p . the pre- 
scribed function a, being independent of the time variable, can be regarded as an 
independent unknown of the following extended version of (|2.5p in the extended 
variable u := (u,a) 

dtu + dxi{u,a) -g{u,a)dxa = 0, 

dta = 0. 

Assuming from now on that the matrix Z?uf (u, a) is diagonalizable for all u G flu 
and a £ R with real eigenvalues Ai(u, a),. . . , A„(u, a) and a full set of eigenvectors 
ri(u,a),..., r„(u, a), it is obvious that (|2.9p admits the same eigenvalues plus 
(with multiplicity 1). Moreover, it admits a full set of eigenvectors if and only 
if Xj{u,a) ^ for all j = 1, . . . , n. In general, (|2.9p is solely weakly hyperbolic 
and, due to possible resonance effects, difficulties arise even in tackling the simplest 
initial value problem, i.e. the Riemann problem; see the pioneering work of Isaacson 
and Temple [3T], as well as Goatin and LeFloch [53] for a general Riemann solver. 
In the rest of this paper, our assumptions will explicitly exclude the resonance effect 
in solutions under consideration. 

While a rigorous definition of nonconservative products will wait until the follow- 
ing section, we can here already provide some preliminary discussion, based on an 
observation by LeFloch [37] for the nozzle flow equations and on the presentation 
in Bouchut [T^] for more general systems. 

With this non-resonance assumption enforced, we then observe from (j2.9p that 
the variable a is a Riemann invariant associated with the eigenvalue A„+i (u, a) := 0. 
In other words, a stays constant across waves associated with any other (non- 
vanishing) eigenvalue and, consequently, the nonconservative product g(u, a)9a;a 
only needs to be defined for {n + l)-contact discontinuities. 
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The entropy inequality (|2.6p should be satisfied as an equality in the sense of 
distributions across standing waves, hence 

(2.10) J"(u+,a+) - J'(u_,a_) =0. 

From the physical viewpoint, we can further investigate the validity of p.lO|) . ob- 
tained as a direct consequence of the augmented form (j2.9l) . To that purpose, we 
specialize (|2.10p first to the case of the shallow water equations (|2.ip 

(2.11) ^ + e(p,) + ^ + a, = + e(p_) + ^ + a_, 
2/3+ P+ 2p± p- 

where m —: p^v^ = denote the mass, and second in the case of the nozzle 

flow equations 



(2-12) TTT-T + <P+) + — = TZ-T-T + <P- 



+P+ P+ 2a%pi_ p_ 

in which m :— a-p-V- = a+p+?;+. The above equations can be implicitly solved 
in p+ away from resonance (see [El EH SSI Hi] for details) and stand at the very 
basis of the design of well-balanced schemes. 

Furthermore, in concrete experiments with fluid flows, for instance in nozzles, it 
is observed that an abrupt change (modeled therefore by a discontinuity) in either 
the topography function, the duct cross-section, or the porosity function, gener- 
ally produces fine-scale features in solutions which may enter in competition with 
complex dissipation phenomena such as friction. To account for such dissipation 
mechanisms, the entropy law (|2.2p or (|2.4p cannot be any longer expressed as a con- 
servation law across the standing wave. The associated entropy dissipation rates are 
the so-called "singular loss of momentum" used by engineers and well-documented 
in the applied literature. It is necessary, on the ground of physical experiments, to 
replace (|2.12l) by the more general condition 

(apv)+ — (apv)^ = 

(2-13) 

J = -(apv)- k(u_, a_), 

with 

J - iapvU (f + eip,) + ^) - (ap.)_ + e(p_) + ^), 

where the prescribed function k : fl^ x R ^> M_|_ defines mathematically the singular 
loss of momentum. The extension relative to (j2.1ip is completely similar. 

We continue this section with a more sophisticated model of compressible flows, 
describing multi-Ruid mixtures, where the variable a introduced previously now 
stands for a fluid mass fraction. Following the celebrated review by Steward and 
Wendroff [54] , (stratified) multi- fluid models may be regarded as two distinct fluids 
evolving with distinct velocities and distinct thermodynamic properties, each prop- 
agating within "nozzles" whose cross sections denoted by a G (0, 1) and 1 — a{x, t), 
respectively depend on the spatial as well as the time variables (see [2] for instance) . 
For notational convenience, a is traditionally denoted by ai (void fraction of the 
fluid 1) and (1 — a) by a2 (void fraction of the fluid 2), with 

ai{x, t) + 0L2{x, t) = 1. 

Furthermore, the evolution of a is now described either via an algebraic closure 
equation (based on the isobaric assumption; see [54] for details) or by considering it 
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as an independent variable governed by a supplementary evolution equation. From 
the point of view of the present paper and in order to avoid instability issues due 
to lack of hyperbolicity of the model, we adopt the second strategy, following here 
Ransom and Hicks 51, and Baer and Nunziato [5] . This approach was extensively 
investigated in recent years; see Gallouet, Herard, and Seguin [22 , Berthon and 
Nkonga jUj, and the references therein. 

In turn, the multi-fluid model under consideration takes the form 
(2.14) 

dtai + Vi{u)dr,ai = \{p2 - pi), 
dt{aipi) + d^iaipi ui) = 0, 
dt{aipi ui) + dx{aipi u\ + aipi) - P/(u) 9^,0:1 = A (u2 - ui) + e9x(Mi d^ui), 

dt{a2P2) + dxia2P2 "2) = 0, 
dt{a2P2U2) + dx{a2P2u\ + a2P2) - Pii^) dxa2 = -A (w2 - ui) + edx{p2 dxU2), 

where u := (ai, aipi, aipiui, a2p2, 02^2^2) is the vector-valued unknown. Here, 
the barotropic pressure laws pi = Pi{pi) are assumed to satisfy the monotonicity 
condition p'^ipi) > 0. The relaxation parameter A > may take arbitrarily large 
values, depending of the multi-fluid flow regime under consideration, while e > 
(the inverse of a Reynolds number) is usually small. Moreover, the (smooth) 
functions Vi : ftu R and Pj : flu R represent the interfacial velocity and 
interfacial pressure, respectively. Following the original work by Ransom and Hicks 
|51j , one can set for instance 

(2.15) Vi{u):=^{u,+U2), Piiu):=^{p,+p2). 

It turns out that, independently of the precise form of the constitutive equations, 
the system (|2.14p admits five real eigenvalues, i.e. 

Vf(u), u,±c,{p,), 

where c'f{pi) := p'{pi) > (as well as a basis of right eigenvectors) if and only if 

(2.16) \Vi{u)-u,\^c,{p,), 1 = 1,2. 

In other words, like the (much simpler) model (|2.9p . the principal (first-order) part 
of (|2.14p is only weakly hyperbolic if (|2.16p is violated. Here again, we tacitly 
assume that solutions under consideration do not develop resonance phenomena. 

One key constraint that arises in choosing the required closure laws for V7(u) 
and P/(u) is the existence of a mathematical entropy pair associated with (|2.14l) . 
Interestingly, the total energy 

U := aipiEi{u) + Q!2/92^^2(u) 

2 

with Ei{u.) :— ^ + ei{pi) \s an entropy for (|2.14p if and only if the interfacial closure 
laws Vi (u) and Pi (u) satisfy the interfacial compatibility condition 



(2.17) 



Vl{u){p2 -pi) +P/(U)(U2 -Ui) =P2Ml ~PlU2 
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for all states under consideration (see |22] and [181 instance) . Indeed, under the 
assumption (|2.17|) . smooth solutions of (|2.14|) satisfy the entropy balance law 

dtU{n) + d^F{vi) = -\{U2 - uif - X{p2~Pif - V, 

U{u) (aipiEiiu) + a2P2£'2(u)) 

(2.18) / N 
J'(u) := UaipiEi{u) + aipi)ui + {a2P2E2{u) + a2P2)u2j , 

2?(u) e^i{dxUi)^ + e/i2(9j;U2)^ - edxipiaid^ui + p.2a2dxU2)- 

The dissipation V formally converges to a non-positive measure when e — and/or 
A +CXD, so that in this limit we do have the entropy inequality 

(2.19) dtU{u)+d.^T{\i) <Q. 

We conclude this section with another setting for complex compressible materials 
which naturally gives rise to hyperbolic equations with viscous perturbations in non 
conservation form. The models under consideration may be regarded as natural 
extensions of the usual Navier-Stokes equations. Such extensions make use of N 
independent internal energies (ei)i<i<Ar for governing N independent pressure laws 
(pi(r, ei))i<j<jv These PDE models take the generic form: 

dtp + da:{pu) = 0, 

N n 

(2.20) dt{pu) + dx(pu^ + y^pi(r, e^)) = e^^^f y^M»(T, ei)dxu], 

i=l 1=1 

dtipei) + dxipuci) +Pi{T, ei)dxU = eMi(T, ei){dxuf , 

where p > denotes the density, u G M is the velocity, and r = l/p>Ois the 
specific volume. Here, e > stands the inverse of the Reynolds number. 

Several models from the Physics actually enter the proposed framework and can 
be distinguished according to the precise definition of the constitutive closure laws 
for the pressures and the viscosities. Precise assumptions on the required state 
laws will be addressed in Section U devoted to the analysis of the traveling wave 
solutions of (ESni)- 

Models from the plasma physics where the temperature of the electron gases 
must be distinguished from the temperature of the other heavy species, typically 
take the form (|2.20p with N = 2 (see [19] for a presentation). Models from the 
physics of compressible turbulent flows can also be seen to fall within the frame of 
PDEs ((2:201) . We refer the reader to [H [HI El [ISl [H] for the mathematical and the 
numerical analysis of several models ranging from two distinct internal energies, 
the so-called laminar and turbulent ones, to > 2 different energies to account 
for a refined description of the turbulent energy cascade. We also emphasize that 
multi-fluid models, as those studied in |11) . enter the proposed framework. 

In most if not all the applications to complex compressible materials, the inverse 
of the Reynolds number e modulating the strenght of the viscous perturbation is 
a small parameter. Solutions of interest therefore exhibit stiff zones of transitions, 
namely viscous shock layers and boundary layers. Viscous shocks cannot be prop- 
erly resolved for mesh refinements of practical interest and we are thus led to study 
the limit e ^ 0+ in the system (|2.20|) . 

There exists several ways to tackle the system (|2.20|) in the limit e ^ 0+ depend- 
ing on suitable change of variables. It can be seen that p.20p can recast equivalently 
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as: 

(2.21) A(w^)5tw^ +^i(w^)a,w^ = e5,(P(w^)5,w^), 
with ^0 regular, or 

(2.22) dtn' + d^T{u') = en{n\d.,n\d.,^u'). 

Namely in (14. 15^ . the diffusive operator writes in conservation form while Aq{w) 
and ^i(w) are not Jacobian matrices of some flux function. By contrast in (|2.22l) . 
the first order operator stands in conservation form but not the regularization terms. 

The precise definitions of the change of unknown w and u is addressed in Sec- 
tion We just highlight at this stage that concerning the equivalent form (I2.2ip 
and provided that suitable estimates on the sequence of solutions We hold true, 
the right-hand side is expected to vanish in the limit e ^ 0+ in the usual sense of 
the distributions. By contrast the left-hand side in non conservation form may be 
handled thanks to the theory of family of paths introduced by LeFloch jS^, Dal 
Maso, LeFloch and Murat ,21^. As far as the next equivalent form ()2.22l is con- 
sidered, the left-hand side now stands in conservation form and can be treated in 
the usual sense of the distributions. In opposition, the right-hand side cannot any 
longer be expected to converge to 0, generally speaking, but merely to a bounded 
Borel measure concentrated on the shocks of the limit solutions. The next section 
provides a convenient framework for handling the required passage to the limit in 
the PDEs ^TF^. 

3. Kinetic relations for nonconservative systems 

Having in mind the examples described in the previous section, we present one 
of the main contributions of the present paper, i.e. the concept of kinetic relations 
for nonconservative systems, which allows us to rigorously define certain noncon- 
servative products arising in the applications. 

Recall that weak solutions to nonconservative systems are defined in the class 
of functions with bounded variation (BV). By standard regularity theorems, such 
functions can be handled essentially as if they were piecewise Lipschitz continuous. 
Throughout the present paper and for simplicity in the presentation, we restrict 
attention to piecewise Lipschitz continuous functions and refer to [21] for details of 
the the DLM theory. 

For simplicity in the presentation, we restrict attention to solutions defined in 
a neighborhood of a constant state in which can be normalized to be the 
origin. We denote by Bsg the ball centered at the origin and of small radius 6a > 
0. Dal Maso, LeFloch and Murat's definition is based on prescribing a family of 
Lipschitz continuous paths 4> — 0(s;ito,Mi) £ Bsg (s G [0,1]), which allows one to 
connect any two points uq,ui in Bs^ for some 6i < Sq. In particular, it is assumed 
that 

(3.1) 0(O;uo,ui) = Wo, (t>{l;uo,ui) ^ ui. 

(See [211135] for the precise conditions, omitted in this short review.) As proposed 
in LeFloch [37) . this family of paths should be determined from traveling wave 
solutions of an augmented model. 

Indeed, it has been recognized that weak solutions u of (jl.l|) depend on the effect 
of small scales that have been neglected at the hyperbolic level of modeling, but 
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are taken into account in the augmented version 

(3.2) dtu + A{u)dxU = R{u,eu^,e^Uj,x, • • • ), 

where i?' = if e = 0. The family of paths determined by travehng wave trajectories 
precisely yields the "missing information" required to set-up the hyperbolic theory. 

A (piecewise Lipschitz continuous) function u = u(x, t) is called a weak solution 
of the nonconservative system if u satisfies the equations (II. ip in a classical 
sense in the regions where it is Lipschitz continuous and, additionally, the following 
generalization of the Rankine-Hugoniot jump relation holds along every curve of 
discontinuity of u. Precisely, for any shock wave connecting two states uq, ui at the 
speed A ~ A(mo,wi), 



uo, X < At, 
Ml, X > At, 



(3.3) u{x,t) 
we impose the generalized jump relation 

(3.4) -A(ui-uo)+/ A{(j){s;uo,ui)) ds(p{s;uo,ui) ds = 0. 



Jo 

Note that, in the conservative case when A{u) = Df{u) for some flux-function /, 
this relation reduces to 

-A (mi - Mo) + /(mi) - /(mo) = 0, 

which is independent of the paths (p and is nothing but the standard jump relation. 

Based on the above definition, one can solve [5T| the Riemann problem for (|l.ip , 
corresponding to the piecewise constant initial data 



(3.5) u(x,0) 



M;, a; < 0, 
Mr, a; > 0, 



where m/ , Ur are constants in Bs^^ with 82 < (5i ■ This construction generalizes Lax's 
standard construction for conservative systems '341 [35] . Recall that (admissible) 
shock waves must be constraint by Lax shock inequalities (for some j — 1, . . . , N) 

(3.6) Aj(mo) > A> Aj(ui). 

The Riemann solver can then be used to design numerical schemes for the approx- 
imation of the general initial value problem, e.g. Glimm or front tracking schemes. 

In certain applications, it has been found convenient to avoid introducing the 
whole family of paths (j). It is precisely our purpose in the present paper to introduce, 
for a particular class of nonconservative systems, a new definition of weak solutions, 
which imposes Rankine-Hugoniot jump relations in the form of "kinetic relations" 
and does not require the knowledge of any "internal structure" for shock waves. 

We will assume that the nonconservative system under consideration formally 
admits N conservation laws, so we consider the system 

(3.7) dtu + d.J{u)^0, 

which consists of conservation laws valid for smooth solutions, only. Our goal is 
to describe singular limits (j3.2p . where i? = i?' is a nonconservative regularization. 
Precisely, we are going to supplement p.7p with N jump relations, referred to as 
"kinetic relation" , which determines the dynamics of shocks in weak solutions to 

(EH). 
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We suppose that in an open and convex domain U C of the phase space, the 
system is strictly hyperbolic, with eigenvalues Ai(m) < . . . < \n{u) and basis 
of eigenvectors U (u) , (u) . Let £ C M be a compact set containing all speeds under 
consideration in the problem. 

Definition 3.1. A kinetic function is a Lipschitz continuous map ^ : U x C ^ M.^ 

satisfying (for j ^1,...,N) 

A, (u)) = 0, u£U, 
\lj{u)-dA^{u,A)\<ci\A-\j{u)\, {u,A)eUxC, 

for some constant ci > 0. Given a a kinetic function <&, a piecewise Lipschitz solu- 
tion u = u{x, t) is called a ^-admissible weak solution to (jl.ip if the differential 
equations (j3.7|) are satisfied in each region of continuity of u and moreover, along 
any curve of discontinuity of u, connecting some values u_,u^ at the speed A, the 
following kinetic relation holds 

(3.9) -K{u+- u^) + f{u+) - fiu^) = A). 

In certain applications, it may be more convenient to express the kinetic func- 
tions in terms of the left- and right-hand states, that is, $ ~ ^{u-,u+). In the 
applications, the kinetic function $ should be determined from traveling wave so- 
lutions of a specific system p.2[) and should be thought of as a "correction" to the 
standard Rankine-Hugoniot relation. 

By introducing the Borel measure denoted by /i* that vanishes in the regions 
of continuity of u and has the mass <I>(m_, A) along its curves of discontinuity, we 
easily see that Definition 13.11 is equivalent to the requirement (see [39] ) 

(3.10) dtu + d.,f{u)^tit, 

which is regarded as an equality between bounded measures. Note that we recover 
the usual conservative case by simply choosing both $ and /i* to vanish identically. 
In general, /i* depends strongly on the function u. 

In the rest of this section we study the case of genuinely nonlinear systems. 
This assumption allows us to use the shock speed as a regular parameter along the 
(generalized) Hugoniot curve. 

Theorem 3.2 (Riemann problem for nonconservative systems with kinetic rela- 
tions). Suppose that (j3.7p is a strictly hyperbolic system in a neighborhood Bsg of 
the origin and admits genuinely nonlinear characteristic fields only, i.e. 

(VAj-r,)(0) >0, j = l,---,7V. 

Let $ = <I>(u, A) be a (Lipschitz continuous) kinetic function defined in the neigh- 
borhood Bs„ X C for some sufficiently small 6 > by 

C C, :-(A,(0)-^,A,(0)+5). 



Then, there exists Si < Sq such that the Riemann problem (|3.5p . (j3.7p with data 
ui,Ur G admits a unique ^-admissible weak solution in the class of piecewise 
smooth solutions consisting of a combination of rarefaction waves and shock waves 
satisfying the kinetic relation. Moreover, the corresponding wave curves are solely 
Lipschitz continuous. 
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Clearly, under the assumptions of the above theorem, p.Sp implies that, for 
j = l,...',7V, 

(3.11) |;j(u)-$(u,A)| <C2|A-Aj(w)p, {\,u)(^UxCj 

for some C2 > 0. 

Proof. Wc want to generalize the proof given in |21j for general families of paths; 
see also the related proof in 28 for nonclassical shocks. We are going to show that 
the given set of jump conditions p.9|) suffices to determine a (generalized) Hugo- 
niot curve uniquely, and we will investigate whether its tangency and regularity 
properties. The rest of the proof (selection of the admissible part of the Hugoniot 
curve, actual construction of the wave curves, Riemann solution) then follows as in 
pi] and will be omitted. 

We denote by Aj(uo,ui) and lj{uo,ui) the eigenvalues and left-eigenvectors of 
the averaged matrix 

A{uo,ui) := / Df{uQ + m{ui—uo))d'm. 



In a neighborhood of the point {uq, Xj{uo)), we consider the kinetic relation 

(3.12) G(A, Ml) := -A {u^ - uo) + .f{u,) - f{uo) - $(uo, A) = 0. 

Fix some index i and let us restrict attention to the (nonlinear) cone-like region K 
determined by the two conditions on ui G Bs^ 

\{ui - Uo) ■ k{uo,ui)\ > C, |A - Xi{uo)\, 
\ui - uol + |A - Xi{uo)\ < (52, 

where a condition on C, > will be imposed below. Observe that G(uo, Ai(Mo)) = 0. 
Multiplying the generalized jump relation (j3.12[) by li{uQ,ui) we find 

= li{uo, ui) ■ {A{ua, Ml) - A) {ui - uo) ~ li{uo, ui) ■ $(uo, A) 

= {\^{uo,Ul) - A)7j(mo,Wi) • {Ui - Uq) -li{uQ,Ui) ■ $(lio, A). 

Therefore, we can express the shock speed A = A(mo, ui) in the form 

(3.13) = A - Ai(uo,ui) + =— — =: il(ui,A). 

li[uo,ui) ■ (ui - Uo) 

Now, observe that the function fl satisfies 

li{ua,ui) ■ (9a$(wo, A) 



^(.i,A)-l 



li{uo, Ul) ■ {ui - Mo) 
> _C2 0(1)'"^~"°''^'^~^'*'"°'" 



\li{uo, Ul) • (ui - Mo)| 

where the constant 0(1) depends only on the fiux. Hence, we have 

g(«i,A)^l + g0(l), 

which is positive provided ci is sufficiently small. As a consequence, the implicit 
function for Lipschitz continuous mappings applies and shows that the implicit 
equation (|3.13p determine the shock speed A = A(uo,ui) uniquely. 
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Next, we consider the remaining components, corresponding to j =/= i: 

lj{uo,ui) ■ $(uo, A) 



uo) 



A{uo,ui) - Xj{uo,Ui) 

Denoting by L{uq) the N x (N — 1) matrix of vectors ^j(wo) for j 7^ h we can 
compute the differential of H, as follows: 



DH 



(mo, Ml) = L{ua) + 0(1) \u 



Mo 



0(1)C] 



|A(uo,ui) - Ai(Mo) 



|A(uo,ui) - Aj(uo) 
0{1)Ci\K{uq,u{)-\{uq)\ 



dA 

^(uo,«i) 



0{l)Ci\A{uo,ui)-X,{uo)\' 



dA 

^(uo,ui) 



2 VA,(wo) 



where we used that A{uq,ui) — Xj{uQ,ui) is bounded away from 0. Hence, we find 



dH 

dui 



(mo,mi) = L(wo) + 0(1) + 0(1) 



aA 

dui 



(uo,"i) 



On the other hand, the ui-derivative of the shock speed satisfies 
OA 



dui 



{uo,ui) 



- VA, wq + O 1 777 7 

2 \li{uo,ui) ■ [ui 



Oil) 



|A-A.(mo)P 



|/i(wo,ui) 
which shows that 



(ui - -uo)|^ 



dA 



-Wo) I 
0(1) 



|A- A^Mo)! 



dA 



\k{uo,ui) ■ (ui - uq)\ dui 



(Uo, ui) 



1 



-(uo,Mi) = - VAi(-uo) + 0(1). 
oui 2 

In conclusion, ^^{uq,ui) — L(uq) + o(l), and the implicit function theorem 
applies to the set of equations H(uo,ui) — 0, which therefore determines a unique 
shock curve s 1— > mi = ui(s;uo), defined locally near uq. Near the base point 
u(0) — Uq, the tangent of this curve is defined almost everywhere and, due to the 
smallness of the constant ci in p.Sp . takes its values in a small neighborhood of the 
eigenvector r.i(uo)- D 

We now introduce a class of nonconservative system to which the framework in 
the previous subsection can be applied. We assume that the first N ~p equations in 
(jl.ip take a conservative form while the remaining p equations are nonconservative. 
In other word, we set u = (w, w) and we consider the nonconservative systems 

dtv + d^g{v,w) = 0, 

dtw + B{v,w) dxV + C{v,w) dxW = 0. 

Here g = g{v, w) G R^"^ while B = B{v, w),C — C{v, w) are p x {N — p) and pxp 
matrix-valued mappings, respectively. 

It must be stressed that the assumption made here refers directly to the set of 
equations listed in or to linear combinations of them. Of course, nonlinear 

functions of the original variable u cannot be considered at this level of the analysis, 
in general, since discontinuous solutions are sought. 
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Our second assumption is the existence of p matlieniatical entropy pairs. That 
is, we assume that there exist k strictly convex functions Uk = Uk{v,w) together 
with their associated flux ~ Fk{v, w) such that 

(3.15) dtUk{v,w)+dMv,w)^0, k = l,---,p, 

holds for all smooth solutions to p.8|) . We search for solutions satisfying the entropy 
inequality 

(3.16) dtUk{v,w)+dMv,w) <0, k = !,■■■, p. 
Many of the models of interest take the form p.8p - (l3.16l) . 

Definition 3.3. Nonlinear hyperbolic systems in nonconservative form that have 
the structure p.l4|) . admit at least p mathematical entropies, and satisfy the non- 
degeneracy condition 

(3.17) det {VM{v, w),-- - , V^,Up{v, w)) ^ 

are called nonconservative systems endowed with a full set of entropies. 

We now focus on the entropy dissipation associated with the entropies Uk- The 
basic idea is to replace the nonconservative equations in the system p.8|) with 
conservative equations for the entropy dissipation but the latter involving a measure 
source-term. Of course it is necessary for {v,w) i— > (v^U{v,w)) to define a change 
of variable, say 11^ j^O- 

Observe first that the inequality ()3.16|) implies a constrain on shock waves, i.e., 
with the notation introduced earlier in p.3|) . 



EkiA;uo,ui) :^ -k{Uk{ui) - Uk{uo)) + Fk{ui) - Fk{uo) 

^^•^^^ <0, 

for all fc = 1, . . .p. On the other hand, the first N —p equations in (2.6) yield N —p 
jump relations in the fully explicit form 

(3.19) - A (ui - vq) + g{vi,wi) - g{vo, wq) = 0. 
Since p jump relations are "missing" , we supply it in the form 

(3.20) £;,(A;uo,ui) = $,(A;uo) < 0, i = l,...p, 

which we refer to as a kinetic relation and where $ is a given "constitutive" function, 
called a "kinetic function" , to be determined case by case in the examples. 

Definition 3.4. Let $ = (0, . . . , 0, $i, . . . $p) be a kinetic function. A piecewise 
Lipschitz continuous function u = [v, w) is called a ^-admissible solution of the 
nonconservative system p.Sp if it satisfies the equations in a classical sense in the 
regions of continuity and if each propagating discontinuity satisfies the N — p jump 
relations (|3.19|) together with the kinetic relations p.20|) . 

We reformulate the main result, in a slightly weaker form which is adapted to 
the present context, since it is natural to assume that the entropy dissipation is of 
cubic order near the base point. 

Corollary 3.5 (Rieniann problem for nonconservative systems endowed with a 
full set of entropies). Consider a nonconservative system endowed with a full set 
of entropies. Suppose that the system is .strictly hyperbolic and genuinely nonlinear 
in the neighborhood of some state u, = (?;*,?«*). Let $i = $i(uo,ui) be a regular 
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function defined in the neighborhood of each speed Xj{u„) for j = I,-- - ,N and 
satisfying for all wq, ui 

^i{uq, Xi{uo)) = 0, 
^^■^^^ 9a$,K,A) = 0(1)(A-A,K))', 

where 0(1) denotes a positive and bounded function. Then, the corresponding Rie- 
mann problem admits a unique admissible solution in the class of piecewise smooth 
solutions consisting of a combination of rarefaction waves and admissible shock 
waves. 

The theory of the present section appHes to the examples listed in Section 2, 
at least when resonance effect is avoided. It is straighforward to include linearly 
degenerate characteristic fields provided the kinetic function is chosen to vanish 
identically for those fields. 



4. Multi-pressure Navier-Stokes system 

In this section, we establish the existence and uniqueness of the traveling wave 
solutions of the multi-pressure Navier-Stokes equations introduced in Section 2, 
under fairly general assumptions on the pressure and viscosity closure laws. The 
equations under consideration were stated in (|2.20l) . Each smooth pressure law 
Pi{T,ei), 1 < i < N , is assumed to obey the second principle of the thermodynamics, 
namely 

(4.1) Ti{T,ei)dsi = dei+ pi{T,ei) dr, 

where Tiir^Ci) > is the corresponding temperature variable and Si > denotes 
the specific entropy. 

The map (t, Sj) i— )> ej(r, s^) is thus well-defined and is assumed to be strictly 
convex. In addition, the following asymptotic conditions are assumed 

(A0\ lim air, Si) ^ +00, lim ei(T, s^) = +oo, lim e,j(r, s,j) = 0. 

\ I T— >0+ Si— >+oo r— ^ + oo 

It follows that 

de de ■ 

(4.3) Pj(t, Si) = — -^(r, Si)>0, Ti{T, Si) = —^{t, s^) > 0. 

OT OSi 

Furthermore, the following assumptions are introduced for any given t > 0: 

dpi 
dsi 

I 

N 

(4.6) lim y^pi(T,s) = -l-oo, lim y^pi(T,s) = 0, 

i=l i=l 

N ^ Nr. 

(4.7) linr V^(r,s) = -oo, lim ^fi(r,s)=0, 

i—1 i—1 

(4.8) lim ^(r,s,) = -^, 



(4.4) |Pl(r,s,)>0, 



(4-5) E|S(-'«)>0, 



1=1 

N N 
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where s = (si, sat). We refer to Menikoff and Plohr [5^ for general properties 
of the fluid equations and the equation of state. Next, the viscosity laws are given 
smooth functions with 

N 

(4.9) fJ.i{T,Si) >0, l<i<N, and ^(r, s) ^ ^,(t, s;) > 0. 

i=l 

To shorten the notation, the PDE's system ()2.20|) is given in the condensed form 

with V in the phase space: 

= {v = (p, pu, {pe,)i<,<N) e p>0,pueR,pe,>0,l<i< N}. 

The basic properties of (j2.20l) are summarized in the following statement. 

Lemma 4.1. The underlying first- order part from h2. 20\) is hyperbolic in fl^ o.nd 
admits three distinct eigenvalues 

(4.10) Ai(v) = u- c(v), A2(v) = ... = AAr+i(v) = u, \N+2iy) = u- c(v), 
(4.11) 

where we set 

(4.12) e^(v)=f;-r*(r,s.). 

1=1 

The extreme fields are genuinely nonlinear while the intermediate ones are linearly 
degenerate. Then, smooth solutions of \2.20]) satisfy the additional conservation 
law 

dtipEY+d, (^{pE}{w^)+Y^^p,{T^,stX^ = ed, (j2p,{T\stXd,u^^ , 
where the total energy reads 

(4.13) (^pE)=^-P^+f2pe.- 

At last, the smooth solutions of \2. 20]) obey the N balance equations 

(4.14) dt{ps,Y + 9. {{ps.Yu^) = e ^'^l[''h d.u^ f. 

As already claimed, changes of variables with distinctive features allow to recast 
(|2.20p either in the equivalent form: 

(4.15) A(w')atw' +^i(w')a,w' = e9,(P(w')a,w'), 
with ^0 regular, or in the form 

(4.16) atU^ + 9,.F(u^) = e7^(u^a,u^a,,u'). 

We briefly discuss the changes of variables involved in (|4.15p and (|4.16[) . Concerning 
(|4.15p . we first observe that summing the N governing equations for the internal 
energies yields 

N 

dtpe + dxpeu + ^^pidxU = ep{T, s){dxu)'^, 
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SO that the foUowing identities are easily checked: 

N 

(4.17) p,i{T,Si){dtpe + dx{peu) + ^Pid^uj = 0, 



N 

l<i< N ~l. 



Since pe = pE — (pu)^/(2p), the conservation laws for p, pu and pE supplemented 
by the {N — 1) balance equations in ()4.17p can be seen to give the equivalent form 
stated in (|4.15p when defining w = {p,p,pE,{pei)i<i<N-i)- A direct calculation 
shows that detAoiw) = ^(r, s)"~^ > 0. 

Concerning system (|4.16p . several changes of variables can be used and we 
advocate in the sequel the change of variables v G fiv ^ u(v) G il^, with 

= {p,pu,{pSi)i<i<N)- 
As underlined in Section [31 both approaches rely on the study of the traveling 
wave solutions of (|2.20p . Due to the frame invariance properties satisfied by the 
PDE model (|2.20p . it suffices to analyze traveling waves solutions associated with 
the first extreme field. With this respect, the main result of this section is as follows. 

Theorem 4.2 (Traveling wave solutions to the multi-pressure Navier-Stokes sys- 
tem). Consider the multi-pressure Navier-Stokes system p.20p when the pressure 
satisfies the positivity, convexity, and asymptotic conditions (|4.2p - (|4.9|) . Let ul G 
riu and a G M 6e given such that 

(4.18) ^^>1' c'{nL) = J2-lTi-L,is.)L). 

c(ul) ^ dr 

Then, there exists a unique traveling wave solution to (|2.20p issuing from the left- 
hand state and reaching some right-hand state U/j G with 

(4.19) 0<^^<1 

The proof of this result will follow from the characterization of a positively 
invariant compact set of 57u- Then the Lasalle invariance principle applied in con- 
nection with a suitable Lyapunov function ensures the existence of a traveling wave. 
Uniqueness is obtained as a simple consequence of the center manifold theorem. 

We gather here some of the notation used repeatedly hereafter and give the pre- 
cise form of the autonomous system which governs the viscous profiles we study for 
existence. Simple but useful geometrical properties induced by the corresponding 
vector field will be then put forward. 

Due to Galilean invariance, it suffices to consider the case of a null velocity a. 
The precise form of the PDE system governing the traveling wave solutions then 
follows when restricting attention to solutions which depend solely on x: 

{pu)x = 0, 

(4.20) (fiu^ + p(t, s))x = (/i(r, s)u^)a:, 
Ti{T,Si){pSiu)x = /ii(r, s)(uj;)^, 1 < i < iV. 

The first equation in (|4.20p implies that the relative mass flux pu has a constant 
value denoted hy m = plUl- As already underlined, we focus ourselves on traveling 
wave solutions associated with the first GNL field; namely we consider m > 0. 
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Observe that the Lax condition (|4.18p expressed expressed for a null velocity a 
reads 

m > plcl- 

Next by integrating once the second equation in (|4.20p . the identity u = mr 
allows one to derive the following {N + l)-dimensional autonomous system: 



(4.21) 



T = ^ (p{T, s) -p(tl,sl) +to^(t - tl)) := / s J^(t,s), 



fl'^{T,s)Ti{T, Si) 



where dots denote differentiation with respect to the rescaled variable x/m that we 
shall refer with little abuse as a time in the sequel. 

This dynamical system is endowed with the following open subset of ]R^+^ which 
will serve as a natural phase space: 

(4.22) n= {uj := (r, s) £ M^+i; r > 0} . 

To shorten the notation, a given function 5* of t and s is simply denoted hereafter 
by «'(w). 

Recall that the total viscosity is assumed to stay strictly positive over ft. 
Then, the regularity assumptions made on all the thermodynamic and viscosity 
mappings ensure that the vector field in (|4.2ip is continuously differentiable 

The unique nonextensible solution of (j4.2ip with initial data ujq in fl is referred as 
to the flow LOQ-t for the times t in the maximal interval of existence {t~ (luq), ^"^(wo)). 
The positive (respectively negative) semi-orbit 7''"(a;o) (resp. 7~(wo)) classically 
denotes the set of states luq ■ [0,i+(a;o)) = {^^^o ■ t : < t < i+(a;o)} (resp. luq ■ 
{t~{ujo),0] — {ojQ ■ t ■ t~{^o) < t < 0}), the orbit being then defined as 7(1^0) = 
7~(a;o) U 'y'^{uJo)- At last, for each ojq in fi, the positive limit set (the so-called 137- 
limit set in what follows) of finds the definition uj^loq) := nt>o7+(wo ■ t), such a 
set is thus empty as soon as t+(a;o) is finite. 

Before we enter the central part of the analysis, let us underline that the (A^-|- 1) 
constitutive variables of (I4.2ip are necessarily kept in their evolution in time in a 
A^-dimensional sub-manifold of i7, the latter being entirely prescribed by the choice 
of the initial data cjq £ i^. This is the matter of the following statement which 
essentially reflects the conservation property met by the total energy ()4.13|) . 

Proposition 4.3. Let loq be a given state in fl. Then the flow ujq ■ t satisfies for 
all time in its maximal interval of existence: 

(4.23) niiOo-t)^n{uo), 
where the regular mapping 7{ : Q ^ W is defined by 

niio) = e(w) - e{ujL) - — (t^ - r|) + (mV^ + p(wl))(t - tl). 

Proof. All the flows under consideration are at least continuously differentiable in 
their maximal interval of existence. The additional conservation law (|4.13p for the 
total energy therefore applies and its differential form reads 

(4.24) {{E + Tp{u))pu]^ = [^^{u^),]^. 
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In view of the algebraic invariant pu — m, ()4.24p once integrated for a prescribed 
wq in VL between time zero and a given time t in (i^ (wo), ^^(^0)) can be seen to 
read 

{E + Tp{i^))pu} it) ~{E + Tp{uj))pu} (0) = Mpt)} (t) - MflT)} (0), 

= {r^H}W-{r^H}(0). 

Since E writes ?77?r^/2 + e{uj), the definition of J- given in (|4.2ip easily yields the 
required identity ()4.23p after some rearrangements in the terms while subtracting 
for convenience to both sides the constant el + m^r|/2 + tlpl- O 

The above statement clearly implies that all the possible heteroclinic orbits of 
(|4.21l) which connect the critical point wl in the past are only made of states w 
such that 

(4.25) -HH = •H(wl) = 0. 

To end up with these preliminary remarks, we point out an elementary but useful 
geometrical property of the flows associated with (|4.2ip which will put restriction 
on possible right connecting states. 

Lemma 4.4. Let ujq be given in 51, then the subset of defined by 

(4.26) n{uJo) ^{uen;s> Sq, H{io) = ^^(wo)} 
is positively invariant. 

The invariance of this region with respect to all positive semi-flows immediately 
follows from the non-negativeness of the A^^'^-last components of the vector fleld 
entering the definition of (|4.21l) . As a consequence, possible heteroclinic orbits 
connecting in the past must entirely lie in 

(4.27) 17(a;L) = {a;el7;s>SL,H(tj)=0}. 

The region (I4.27|) will play a central role in the derivation of positively invariant 
compact sets. 

Here, we exhibit some important features of the linearization LX(lJc) of the 
vector field X at equilibrium points lJc, i-C- at states satisfying J-{ll)c) — 0. We 
check in particular that such states are always non-hyperbolic points for which the 
space M^+^ writes as a direct sum of the eigenspaces associated with LX{u}c) under 
the following non-degeneracy condition: 

(4.28) drJ^iUc) = + drP{Tc, Sc) ^ 0. 

In that aim, let us state some basic facts concerning the linearization LX{oji^. 
The requirement J^{ujc) = is easily seen to enforce all the partial derivatives of 
the A''*''-last components of X to be identically zero (since these components are all 
proportional to J"'^). Under the nondegeneracy condition (|4.28l) . there consequently 
exists only one non trivial eigenvalue namely drJ^{(^c)/fJ-{'-^c) while A = is a 
semisimple eigenvalue of LX((jJc) of multiplicity N . Furthermore, the corresponding 
eigenspaces T(ujc) and T'^iujc) are respectively the span of ei and 62,..., bat where 
{e-i]i<i<N+i stands for the canonical orthonormal basis of M^"*"^. 

Equipped with these results, the center manifold theorem ensures the existence 
of two locally invariant manifolds >V(a;c) and W"^(a;c) (the so-called center manifold) 
of class at least and C° respectively which go through Wc and are respectively 
tangent to T{ujc) and T^(wc) at this point. The regularity properties above are 
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indeed inherited from the continuous differentiabihty of the vector field, according 
to this theorem. 

Assuming drJ'i'^c) to be positive (respectively negative), i.e. assuming the cor- 
responding sign for the unique nontrivial eigenvalue of LX{ujc)] W(aJc) is classically 
referred as to the unstable (resp. stable) manifold with superscript u (resp. s). Re- 
call that the unstable manifold of a point is the manifold composed of the totality 
of the orbits which tend exponentially fast to the point in negative time; the stable 
manifold being defined conversely. Then by well-known topological considerations, 
two rest points namely ujl and ljij are connected by a heteroclinic orbit 7 precisely 
ii-i (ZW{uJL)r\W{uB). 

An obvious requirement for the existence of a heteroclinic orbit connecting ujl 
in the past is then 

drT{uJL) ^ nri^ + drp{uJL) < 0, 
but the validity of such an inequality is precisely the matter of the Lax condition 
(14.181) . Conversely, a possible connecting point uji in the future is necessarily 
subject to the condition drJ^{ujB,) > 0. 

Now and since the unstable manifold yV"(cji) is one dimensional, there exists 
locally exactly two solutions of (|4.2ip which approach cj^ as i —00. Arguing 
about the property of W"(wl) to be tangent to ei, the associated almost horizontal 
orbits approach lo^ from the two opposite directions t > and t < t^. With clear 
notation, 7>(aji) (respectively j^{ujl)) will denote the first (resp. the second) orbit. 

The following assertion discards the solution converging to ujl for negative times 
from the region t > t^. Note that such a result precisely precluded expansion 
shocks to admit viscous profiles. 

Proposition 4.5. There is no heteroclinic orbit of the dynamical system \4.21\j in 
the domain J\f := {uj ^ VL : t > tl, s > s^}. 

Consequently, only the second solution can give rise to a heteroclinic orbit. Since 
the vector field X : Q, ^ M.^^^ is Lipschitz-continuous, the uniqueness part of the 
celebrated Picard-Lindelof theorem readily gives the following. 

Corollary 4.6. There exists at most one heteroclinic orbit of the dynamical system 
{4.21^ which connects lul in the past. 

The proof of Proposition 14.51 will follow from the following statement. 

Lemma 4.7. Any given state tu distinct from ujl in the set {r > Ti,s > s^}, obeys 
T(uj) > 0. 

Proof. Observe that the positiveness assumption (|4.5|) on the Griineisen numbers 
implies that for all u! in the region under interest (s > s^): 

(4.29) f{T):=T{T,SL)<T{Lu) 

with equality iff s = s^. In particular, for t = tl we have T^t^^s) > as soon as 
s > sl. Next, considering t > tl, the following identity 

f'{T)^m' +drP{T,SL) 

clearly yields, under the assumption (|4.4|) of positive fundamental derivatives, that 
d^^p{T, s) > for all w G rj andn therefore, 

/'(r) > /'(tl) = - {pLCLf > 
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thanks to the Lax condition (|4.18p . It immediately fohows that /(r) > /(tl) = 
J-{ujl) = as soon as r > tl with equahty to zero iff r = TL. The inequality (14291) 
then gives the required conclusion. □ 

Proof of Proposition \4-5\ Let be the unit inward normal at the following hyper- 
surfaces {r = tl,s > sl} and {r > tl,s = s^} for all states in these sets. Note 
that these sets are the lower part of the boundary of Af. In view of the definition 
of the vector field, Lemma [4.71 implies that, for such states, X{uj) • n,^ > 0. As a 
(well-known) consequence, Af stays invariant for all positive semi-flows. The re- 
quired conclusion follows again from Lemma 14.71 which says that no critical point 
exists in Af. □ 

We now stress another important consequence of the local properties of the phase 
portrait at the rest point wl. By opposition to the states in the orbit 7>(aJL); 
the second orbit 7<(a;L) emanating from the region r < is by definition made 
of states to that at least when close enough to but distinct from lul give rise to 
a compression: namely locally J^{uj) < in view of the governing equation for 
T. In turn, this simple observation implies that the viscous profile under study 
must remain uniformly compressive. This claim is a consequence of the following 
statement. 

Lemma 4.8. The following set 

(4.30) I^{ujEn:T<TL, s > SL, ^"(0;) < 0} 

is positively invariant under the action of the dynamical system ^.21^ . 

Proof. The above assertion is trivial for states G X which satisfy J^(wo) — 0. 
Considering states loq with the property J-{ujq) < 0, we observe that the positive 
semi-flow loq ■ t necessarily satisfies J-{ujq ■ t) < for all time in [O,t+(a;o)). Indeed 
assuming the existence of a finite time t^ in this interval with the property J-{ujo ■ 
tc) = would result in a critical point ujq ■ tc for the dynamical system (j4.2ip . But 
by the Lipschitz-continuity property of the vector field in il, it is well-known that 
such a point cannot be reached in finite time. □ 

The orbit 7<(wl) is therefore trapped in the region I. We now establish that, 
in addition, this orbit must remain within a compact subset K of I. This will 
imply that 7<(cjl) is relatively compact. Well-known considerations imply that 
the associated nj-limit set is nonempty, compact and connected. The existence of 
K primary stems from the following result. 

Lemma 4.9. Let loq be a given state in fl. Then the positive semi-orbit ^'^(ujq) 
has no limit point in the set {t — 0}. 

This assertion immediately gives that the orbit ^^(uj^) has the same property. 

Proof. For all time t in [0,i^(aJo)), the positive semi-flow uja ■ t is known to obey 
Hiuj) — HiuJa) < 00 and s > sq. Arguing about the positivity of all the tempera- 
tures Ti = dsfH, we immediately get 

h{T) :=7^(t,So) <H(r,s), 

with the property that /i(t) goes to infinity as r goes to zero (see the asymptotic 
condition (|4.6p ). 
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Assume that 7+(a;o) n{r = 0} is nonempty. As a consequence, for all e > there 
exists e (0,t+(wo)) such that < T|^(,.f^ < e. Necessarily there exists Eq > so 
that /i(T|^Q.t^^) > Ti^uJo) and this rises the contradiction with the preservation of 
T-L{i^q) along the orbit. □ 

We have proven that any given positive semi-flow of the dynamical system (j4.2ip 
with initial data in Vl satisfies t"'"(u;o) = oo (since s > Sq). We now conclude 
with the existence (and therefore uniqueness) of the required viscous profile. 

Proposition 4.10. There exists a state cuji in H^^(O) H J^^^(O) which is connected 
by 7<(wl) in the future. 

Proof. We first establish that the specific entropy vector s stays upper-bounded 
along all positive semi-flows with initial data in I. For flxed r in (0,Ti), the 
conditions (|4.3p and (|4.5p shows that H(t,s) rises arbitrarily with s. The same 
steps as the ones involved in the previous proof, apply to give the required result. 
As a consequence, 7<(a;L) must be included in a compact subset, namely K, of the 
positively invariant region I. This orbit is, therefore, relatively compact and its 
•nj-limit set is non-empty. This limit set must be included in H~^(0). To conclude, 
observe that r, when understood as a mapping of the variable w, trivially yields a 
Lyapunov function on I where, by construction, -F(w) < 0. The LaSalle invariance 
principle applied in connection with this Lyapunov function then ensures that the 
non-empty tu- limit set is included in {w G X : J- (to) = 0}. This establishes the 
existence of ujr. □ 



5. End states for viscous layers with varying viscosity 

The existence (and uniqueness up to translation) of traveling wave solutions to 
the multi-pressure Navier-Stokes equations was established in the previous section 
for N viscosity laws satisfying the non degeneracy condition (j4.9l) . Being given a 
fixed state e and a velocity a according to the condition (|4.18p . we aim here 
at characterizing the subset of D, made of all the states tdR that can be reached in 
the future by a traveling wave with speed a and connecting uj^ in the past. We 
naturally expect the exit state ujn to depend on the specific form of the iV-uple 
of viscosity laws. The dynamical system (|4.2ip shows that such a dependence is 
in the ratios of the viscosity laws. As a consequence, possible states ujji to be 
reached in the future from a fixed ljl (at some speed a) generically depend on 

— 1 degrees of freedom. The set of exit states we are seeking is thus expected to 
have {N — l)-dimension. 

It will be convenient to study the projection of this set onto the following positive 
cone of (understod as the space of the specific entropies s = (si, sn)) with 
origin s^,: 

(5.1) S+{s) = {s e / s St -t- Aa, a e 5^ , A > 0} . 

Here, denotes the (positive) part of the unit sphere in defined by 

Sf = {^eR^; ||a||=l}. 

For all possible entropies s^^ in the cone (|5.ip . the existence simply comes from 
the property that the heterocline solutions of Theorem 14.21 obey > st and a 
strict inequality holds for (at least) one specific entropy Si {1 < i < N). 
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We show hereafter that the projection in the half cone (|5.ip of the states ujr that 
can be reached when varying the definition of the N viscosity laws, is a smooth 
manifold with co-dimension one: 

(5.2) C = {s e 5+(s) / s = Si + Ao(a)a, a £ 5f }, 

for some suitable mapping Ao(a) : a e 5^ h-^ Ao(a) G K which precise definition 
will be given latter on. The derivation of the proposed manifold is performed in 
two steps. In a first step, we analyse closely all the critical points (tc, Sc) of the 
dynamical system (I4.21[) . i.e. the solutions of 

(5.3) .F(rc,sc) = 0, 

without reference to a precise A'^-uple of viscosity laws. 

We emphasize that eligible critical points that can be reached from the state ujl 
in the past must preserve the total energy as stated in (|4.25p . Such states must 
therefore solve in addition 

(5.4) H(tc,Sc) = 0, with Sfl>Si. 

Analyzing the solution of (|5.3|) - (l5.4p will give birth to the manifold (|5.2p . 

In a second step, we will establish that any given value s in the proposed manifold 
can be actually achieved for at least one suitable iV-uple of viscosity laws. As a 
consequence, the manifold (|5.2[) is entirely made of all the specific entropy s^^ that 
can be reached in the future by a traveling wave solution with speed a and issued 
from wl, when varying the definition of the N viscosity laws. 

Let us outline the content of this section. We first analyze the mappings s G 
5+(s) Tjr{s) G R+ that solve 

(5.5) ^(t^(s),s) = 0. 

We then characterize the mapping s G 5'+(s) ^ th{^) £ R+ solving 

(5.6) H(r„(s),s) =0. 

Equipped with these two families of functions we will study for existence values 
of the specific entropy in ^^(s) that satisfy Tjr(Sc) = r-H(Sc), namely values of 
s that simultaneously solve (15. 5p and (|5.6p . We now state the main result of this 
section. 

Theorem 5.1. Assume that |^.g[ )-( P75[ ) on the thermodynamics are satisfied. Then 
there exists a unique map T G C"(/C,M!i_) nC"'^(/C,R^) where K. C <S'+(s) reads 

(5.7) /C = {sg5+(s) /s==SL + Aa, a G 5^ , A g]0, Ao(a)[}, 

for some smooth application Ao G (5:^, M!i_) with the following properties: 

(i) H{T(s),s) =0, for alls eiC, 

(ii) J-"(T(s),s) = 0, for all s ^ C where 

(5.8) C = {sgS'+(s) / s = SL+Ao(a)a, aG5^}. 
In addition T obeys 

(iii) ^(r(sL),si) = 0, 

(iv) J'(r(s),s) < 0, for all s G /C. 
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The mapping Aq : — > M.*^ will be built in the course of the proof. Rephrasing 
the above result, the function T(s) with s £ C, simultaneously makes vanish T and 
H, so that all the values s in the smooth manifold C are candidate for being reached 
in the future from the state cj^ via a traveling wave with speed cr for suitable choice 
of the N viscosity laws. 

We now show that all the values s in the manifold C, defined by (|5.8|) . are actually 
eligible candidates for entering the definition of the specific entropy in exit states 

UJR- 

Lemma 5.2. A state ujl being given in and a velocity cr being prescribed according 
to \4.18^ . For any given s G C, there exists at least one relevant definition of the 
N-uple of viscosity laws which yields a traveling wave solution with speed a issued 
from lul and connecting a state ioj^ in the future with S/j = s. 

Proof. The proof of this result makes use of particular viscosity laws under the form 

(5.9) M,(r,s,) = M?r,(T,s,), A'°>0, l<i<N. 
The non degeneracy condition (|4.9p is satisfied as soon as 

N 

(5.10) E^^>0' 

since each of the temperature law Ti{T, s^) is assumed to be positive. Without lost 
of generality, we assume fi% > 0. 

We stress that viscosity laws which linearly depend on the temperature naturally 
arise in the kinetic theory for dilute gases. We refer the reader to the book by 
Hirschefelder, Curtiss and Bird |29| . 

Observe that viscosity laws in the special form (|5.9p let evolve each specific 
entropy according to 

fi^[T, S) 

We thus infer the following {N — 1) balance equations linking the evolution of the 
first {N — 1) specific entropies Si to the last one: 

Si^l^SN, l<i<N-l. 
Mat 

Since the ratios n^/fi^ are constant real numbers, we deduce: 
,0 \ ,,0 



^sn) iO = 4 - ^sj(,, for all ^e . 



l^N / l^N 

We therefore end up with [N — 1) jump relations 

(5.11) ,f _,i^4(4-4), 1<*<A^. 

We emphasize at this stage that — > in view of our assumption /i^ > 0. 
From the jump relation (|5.1ip . we therefore get 



(5.12) s«-Si = ^i^4^ 

Ma 



V A^°a 



WHY MANY THEORIES OF SHOCK WAVES ARE NECESSARY 



25 



which is obviously in the form Ao(a)a, for a in , given by 

('^")i<,<Ar- 



Next and up to some relabelHng in the viscosity in order to allow /x^^r to vanish, 
any given a G gives rise to an admissible A^-uple of viscosity coefficients. This 
concludes the proof. □ 

Remark 5.3. The identity i5.18\) shows in addition that the mapping a i— > Ao(a) 
can be built as soon as the jump in the last specific entropy — s^ is known. This 
evaluation can be performed numerically; see, for instance, 8 . 

We now give a proof of the main result of this section, namely Theorem 15.11 In 
that aim, and as already claimed, we propose to first study for existence in 
the roots t(s) of T{t, s) = 0. Then, we shall study their distinctive properties by 
investigating the values of ?^(t(s),s). 

Proposition 5.4. There exists two maps, we denote by belonging to C^(T>D 
{sl},RX) nC°{D,R*^), where V is the subset of S+{sl) defined by 

(5.13) V^{seS+{sL) /s==SL + Aa, ae5^, Ae]0,A(a)[}, 
for some A e 0^(5^ ,M+) with A(a) > Ao(a) for all a e , so that 

F{t^{s),s) ^0, set). 
In addition, these two families of roots are interlaced according to 

(i) T-{s) < T+(s) < TL, for all s G V\{sl}, 

(ii) T+(s) = TL in V ifjs = sl with t_^{sl) < t+(sl) = tl, 

(iii) T-(s) = r+(s) inViffs = + A(a)a, a G 5^. 

Again, the map A : S!^ — > M*^ will be built in the course of the proof. But from 
now on, notice that K. C T>. We shall show that for fixed s E V, F{t, s) = only 
admits t='=(s) as roots and cannot be solved in r for values of s in S^{sif)\D. As 
a consequence, all the critical points (r(s),s) of ()4.2ip are necessarily achieved for 
s G P so that r(s) must coincide with either t~(s) or t+(s) for suitable values of 
s G 2? : i.e. such that 7^(t~(s),s) = or H(t+(s),s) = 0. In this way, let us state 
some properties of % with respect to the above two families of roots. 

Proposition 5.5. Using the notation in Propositions \5.i\ and \57^ we have 

(i) H(r+(SL),SL) =0,, 

(ii) ■H(r+(s),s) > 0, for all s G V\{sl}, 
while 

(iii) n{T- {s),s) < 0, for all s G /CU{sl}, 

(iv) -H(t-(s),s) = 0, for all s £ C = {s e V / s = sl + Ao(a)a, a G 5f }, 

(v) -H(t-(s),s) > 0, for all s G 

Put in other words, the critical points of the differential system (I4.2ip necessarily 
coincide with the set {t+(sl),sl} and {((t~(s),s) / s G C}. Keeping this in 
mind, we next analyze the roots r(s) of H(t(s), s). The following claim states that 
Ti admits three distinct branches of roots in IC. A particular attention is paid to 
single out a branch 7" obeying the requirements: 

(5.14) T{sl)^t+{sl) together with T{s) ^ (s) for aU s G C, 
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as put forward in Proposition 15.51 

Proposition 5.6. There exist three maps in C'^{K,R*^_) fl C""^ (/C, K*^) respectively 
denoted by 'T, 7", T" : ^ — > K^, so that: 

(i) n{T{s), s) = H(f (s), s) H{f{s), s) = 0, for all s e £. 
These are interlaced with the roots (s) of J- as follows: 

(ii) f (s) < T-(s) < r(s) < T+(s) < r(s), for all s e IC, 

(iii) r(si) < r(sL) = r+(si) - r(sL), 

(iv) r(s) = r-(s) = r(s) < r(s),, /or all see. 

Observe tlrat the intermediate mapping T fulfills the requirements (|5.f 41) so that 
Theorem 15.11 is established. 

We now give the proofs of Propositions 15.41 to 15.61 Proposition 15.41 relies on the 
following two technical lemma. 

Lemma 5.7. For all fixed s G S^{sl), J-{.,s) admits a unique minimum in r we 
denote f(s) where f G C"'^ (5*^(8^), R!^) with f(s) < T]^ for all s G 2?. This minimum 
obeys: 

(i) J"(f (s), s) < 0, for all seVU {sl}, 

(ii) J"(f(s), s) = 0, for allseT := {s G_5+(sl) / s = sl + A(a)a, a G S^}, 

(iii) J'(f(s),s) > 0, /or all s G S'+(sl)\P, 
where the set T> has been defined in Proposition \5.4\ 



Lemma 5.8. For all fixed s G S^{sl), J-{.,s) is strictly decreasing (respectively 
strictly increasing) for all r G (0,f(s)) (resp. for all t > f{s)) and achieves the 
following limits 

lim J^(r, s) — +00, lim T{t, s) = +oo. 

As a consequence, J-{t,s) — can be solved in r only when s G 2?, with exactly one 
solution when s G F and exactly two solutions t^(s) for s G I?\r. These solutions 
define two maps T± G C°(2?,R!;) nCi(2?U {sl},R!;) with the following properties: 

(1) T"(s) < f(s) < T+(s) < tl, /or s G V, 

(2) T"(s) = f(s) = r+(s) < TL, /or all s = + A(a)a, a G , 

(3) T"(sl) < f(SL) < r+(sL) TL. 



We now establish Lemma 15.71 underlining that the set T) entering the Proposition 
15.51 will be explicitly derived in the course of the proof. 

Proof of Lemma \571[ Let s be fixed in ^^(s). Arguing about the smoothness of the 
internal energies, the map r H> J^(., s) is at least of class C^(R!^). Easy calculations 
then yield for all r > 0: 

dT , , dp , . 2 d"^^ , N d'^P , X 
On the one hand, the map r ^(''"i ^) is strictly increasing in view of the genuine 



nonlinearity assumption (14.41) for the total pressure. On the other hand, assump- 
tions (|4.7p on the asymptotic behaviour of ^ imply that: 

lim — — (r, s) = — oo and lim — — (r, s) = rn^ > 0. 
T->-o+ OT T^+oo dr 
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As a consequence, for all s G S'+(sl), there exists a unique f(s) > so that 
^(f(s),s) = 0. This defines a map f in (S'~'"(s/,), R^) thanks to the implicit 
function theorem. Note that the assumption (|4.18p on the relative Mach number 
implies that drT{TL,SL) > while drJ-{f{sL),SL) = 0, therefore hypothesis (|4.5p 
ensures: 

(5.15) f(si) < TL. 

Next, we construct the set V C 5+(sl) introduced in Proposition l5.5l when studying 
for existence the zeros of s e 5+(si) i-> J^(f(s), s). In this way, we first notice that 
by definition of t(s), for all s G S^{sl) we have: 

T{f{s),s) =p(f(s),s) ~p{tl,sl) - — (f(s),s)(f(s) - tl). 

Introducing the auxiliary function (j) : R+ x — > 5+(sl) defined by 

a) = T{f{sL + Aa), Si + Aa) 
straightforward calculations give 

(5.16) #(A,a)- ^ (^P^^ ' ( ^ ^'P^^^ ' ^"P^ 




dX ' , \ dsi \ St^ dsi drdsi 

l<j<N 

But differentiating the identity (s), s) = — m^, valid for all s G S^{sl), 

easily implies that ()5.16p reduces to 

(5.17) — (A,a)= ^ ^ (f(sL + Aa),SL + Aa) a^. 

l<i<N * 

This derivative is therefore strictly positive for all a G Sf, as follows from (|4.5I) . 
Next, arguing about the strict convexity in r of the total pressure and the property 
(|5.15p expressing that f (sl) < tl, we get 

(j){0, a) = p(f(sL), sl) - p{tl,sl) ~ ^(^(^l), sl)(^(sl) - tl) < 0. 

To conclude, we need to check that for all a G Sf, the map A i— > (f){X, a) achieves 
positive values for finite values of A. Indeed, the implicit function theorem will thus 
ensure the existence of a map A G (5^, R;) well defined in V with the following 
property: 

(5.18) 0(A(a), a) = 0, for aU a G 5^ 

with (/)(A,a) < (respectively > 0) for all A < A(a) (resp. A > A(a)). This is 
nothing but the required result. To establish the validity of (|5.18l) . we show that 
for all a G there exists A* (a) > so that 

(5.19) f(si + A4a)) =Ti. 
Indeed for such values of A, boils down to 

0(A*(a),a) = p{tl,s) - p{tl,sl) > 0, 

as follows from (|4.5p . To derive (|5.19p . we introduce the auxiliary smooth function 
6* G Ci(R+,R) defined for aU a G 5f by 
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For any given a in , we establish the existence of solutions to d{X) — 0, A*(a) 
being chosen to be for instance the smallest one for definiteness. Existence of such 
solution(s) readily follows from the assumption (|4.18l) on the relative Mach number 
ensuring that 9{0) > while the asymptotic condition (|4.7p ensures d{X) < for 
large enough A. Note that the solutions under consideration are strictly positive. 
In addition and since A(a) < A^(a) for all a G , we obtain 

(5.20) f(s) <TL, s e V, 

which concludes the proof of Lemma 15.71 □ 

Proof of Lemma \5.8[ Let be given s in ^^(si). By definition of t(s), J-{.,s) 
achieves the monotonicity properties stated in Lemma 15.81 the required limits im- 
mediately follows from the asymptotic conditions (|4.6|) . The study of the sign of 
J^(f(s),s) (as we proposed earlier) obviously implies that for fixed s e S'+(sl), the 
equation -F(t, s) = has exactly two solutions t~ , r"*" in 'D\T so that < f(s) < 
T"*"; this equation has exactly one solution, namely t(s), when s G F and has no 
solution whenever s e S^{sl)\'D. In addition, using the notation introduced in the 
proof of Lemma ISTtI it can be easily seen that the following limits hold true: 

(5.21) lim r±(sL + Aa) =f(sL+A(a)a), for all a e 

These observations allow for the definition of two maps : V ^ satisfying: 
J'(r±(s),s) = 0, for all s G P, 

and so that 

(5.22) T-(s) < f(s) < T+(s), seP\F; (s) = f(s) = t+(s), s G F. 

Then the above inequalities yield drJ-{T^{s), s) 7^ for all s G P U {s^} in view 
of (j5.22|) so that the implicit function theorem ensures that T± G C^{VU{sl},R*+) 
while (f53T|) gives that t± g C°{V,Rl) . 

Next, focusing to some given s G 'D\{sl}, we observe that 

J"(tl, s) = p{tl,s) - p{tl,sl) > 

so that necessarily, cither tl < t^(s) or t+(s) < tl- In addition, the identity 
F{tl,sl) = expresses that either t^{sl) = tl or T+(si) = r^. But Lemma [5?7l 
ensures that f (s) < tl for all s G P. This concludes the proof. □ 

Equipped with these two lemmas, the proof of Proposition 15.41 is essentially 
completed: the required inequality A(a) > Ao(a) for all a G will be deduced 
from the derivation of the set K. we propose hereafter. We will need the following 
technical result. 

Lemma 5.9. For all s G F, 'H{f{s),s) > 0. 

This statement actually indicates that there is no critical point on F. 
Proof. To shorten the notation, let us introduce 

N 

e('r>s) = ^e,(T,Si), 

i=l 
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and consider the auxiliary function tp £ C^(r,K) defined for all s G F by 

2 

ID 

m = <t(s), s) - e(ri, s) + iJi(f(s) - r^) - — (f (s)^ - r|). 
Here we set 

Hl = rn^TL +p{tl,sl), 

so that H(t, s) recasts as 

H(f(s),s) = VCs) + e(TL,s) - e(rL,SL). 

Arguing about the identity 7^(f(s), s) = valid for all s G F (see Lemma lSTTl fii)). 
we have 

Hl = p{t{s), s) + m^f(s), m2(f(s) - tl) = p{tl,sl) - p(f(s), s), 
which gives successively for all s e F: 

■0(8) = ^{t{s),s) - e{TL,s) + (f(s) - Tl) ( p(f(s),s) + — (t'(s) - tl) 

^ e(f(s),s) - e{TL,s) + (f(s) - tl) ^p(f(s),s) - i(p(f(s), s) - p(tl, sl)) 

Moreover, the two identities J-{f{s), s) = and ^(t, s) = valid for all s G F are 
easily seen to give for the s under consideration: 

dp 

p(t(s),s) -p(tl,sl) = (r(s) - tl) — (t(s), s). 
Consequently, for all s G F 

V'(s) = e(f (s), s) - e(rt, s) + (f(s) - rL)p(f(s), s) - -(f(s) - TL)2^(f(s), s). 
To conclude, we show that 

(5.23) n{f{s), s) = 0(f(s), s) + e(Ti, s) - e(TL, s^) > 0. 

Since ff^(TL,Si) = Ti{TL,Si) > then e(TL,s) — e(rL,SL) > for all s G F since 
sl ^ F. Indeed, observe that Lemma [5.81 implies that equality to zero holds iff 
s = sl but Sl ^ F. 

To show ()5.23p . we study the following auxiliary function ^> € {^+, , setting 
for fixed s G F: 

1 dp 

*W = e(r,s) - €{tl,s) + {t- tl)p{t,s) - -(t - tl)^ — (t,s). 
Easy calculations give 

with ^'(tl) = 0. Consequently, ^(r) > for all r < r^. Since f(s) < tl (see 
Lemma 15. 7p then ^(f (s)) > for all s G F, and we thus obtain the required 
inequality: H(f(s),s) > 0. □ 



Proof of Proposition \5.5\ We first establish the required properties of T-L related 
to the branch of solutions r+. Arguing about the identity t+(s) = f(s) for all 
s G F, the technical Lemma [5.91 allows to restrict ourselves to s G T)\r where t'^ 
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is continuously difFerentiable. For such s, the identity J^(t+(s),s) = re-expresses 
equivalently: 

(5.24) m2(r+(s)-Ti) = (p(rL, s^) - p(t+(s), s)) , . 

Let us evaluate 'H(r+(s),s) as follows: 

H (r+ (s) , s) = e(r+ (s) , s) - e(ri , SL ) + (r+ (s) - tJ p{t+ (s) , s) + ^ (t+ (s) - ) 



-W(t+(s),s) = e(r+(s),s) - e{TL,SL) ~ ^ {p^T+ {s),s) ~p^tl,sl)) 



where e(T, s) = X^il^i ^ii^^ ^i)- Using (|5.24p . we then obtain 

1 

Let us introduce the auxiliary function Q : K.*^ x 2? — > K by setting 

Q{t, s) = e(T, s) - ^^P^{t, s), 

so that for all s G 2?: 

(5.25) H(t+(s), s) = e(r+(s), s) - e(rL, s^), 

with s ^ e(T+(s),s) e Ci(P\r,]R). Since t+{sl) = tl by Lemmas [EH and [52S1 
reads equivalently: 

n{T+{s),s) = e(T+(s),s) - e(T+(si),si). 

Moreover, for all s G 2?\r we have: 

^^(r+(s),s) = — (s)-(.+ (s),s) + -(r+(s),s), 

3t f)f ■ 

= _ (s)^(r+(s),s) + -^(r+(s),s) 

de 

(5.26) = ^(^+(s),s)=r,(r+(s),s)>0, 

where we used the identity 7^(r+(s), s) = 0. Consequently, we deduce that 

e{T+{s), s) - 0{t+{sl), sl) > 0, Vs e p\r 

with equality to zero iff s = (see Lemma 15. 8p . Combining the previous steps 
with Lemma [5.91 gives the required properties (i) and (ii). 

We now derive the remaining properties of Ti related to t^. Observe that the 
technical Lemma 15.91 immediately gives 

(5.27) H(r-(s),s) > 0, s e F, 

since t^(s) = f(s) for the s under consideration. We can now obtain the following 
estimate 

H{t-{sl),sl) < 0. 

To that purpose, let us introduce the following auxiliary function ^ : — > R 
setting: 

^{T)=n{T,SL). 

Since iP'{t) — — 7^(r, s^) for all r > 0, Lemma [5.81 is easily seen to imply that ip 
strictly increases in ^r^(sL), t+(sl)) with H(t+(sz,), s^) = as just established. 
This yields inequality ([S]). 
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To conclude the proof, we follow exactly the same steps as those developed in 
the proof of Lemma [5.71 devoted to the derivation of the subset K. G 5+(sl). 
We introduce the following auxiliary function defined by 

$(A,a) =?^(T-(sL + Aa),SL + Aa), a G 5^^, Ae[0,A(a)[. 

Note that this function is continuously differentiable on its domain of definition 
since, in view of Lemma [5.81 the function (r, A) i— >■ t~{sl + Aa) is differentiable. 
Straightforward calculations then give 



(5.28) 



^(A,a)=-J'(r (sl + Aa),SL + Aa) f V ^— a, 
+ ^ Ti (t"(sl + Aa),SL + Aa) Oi. 



l<i<N 

But the following identity J-{t^ (sl + Aa), + Aa) = holds true by definition for 
all a e 5^ and A £ [0, A(a)] so that reduces to 

— (A,a)= 2^ T,{t (sl + Aa),SL +Aa)aj > 0. 

l<i<A' 

Arguing about the inequalities (|5.27p and ([5|) , the implicit function theorem implies 
the existence of a map Aq G C^(5:5^, M!i_) with the following properties: 

$(Ao(a),a) =0, a G 5^^, 

together with $(A,a) < for all A G [0,Ao(a)[ and $(A,a) > for ah A G 
]Ao(a), A(a)]. This concludes the proof of Proposition 15.51 □ 

Proof of Provosition \5. 6\ Arguing about the identity 9r'H(T, s) — —J-{t,s) valid 
for all (t, s) G x /C, Lemma 15.81 immediately implies that the smooth map 
T n- 7{(t, s), s being fixed in JC, strictly decreases in ]0, t~(s)[ and ]r+(s), +oo[ while 
it strictly increases in ]t~ (s), r+(s)[ with the following limits limT.^o+ T^i^j s) — +oo 
and linv^-oo H(r, s) — — oo in view of the asymptotic conditions l4.2l In addition, for 
all s G /C, we infer from Proposition 15.51 that ■H(r^(s),s) < and ■H(r+(s),s) > 0. 
These observations allow the definition of three maps, namely T, T, T ; /C — > 
with the following properties: 

H(r(s), s) = H(r(s), s) = •H(r(s), s) = o, s g /c, 

together with 

< r(s) < T-(s) < r(s) < T+(s) < f (s), s G /C. 



Next, using the notation introduced in the proof of Lemma 15. 7) we first compute 
for all a G 5f : 

lim T(sl + Aa) < lim T(sl + Aa) = lim T(sl + Aa) = t^(sl), 
\^o+ A^0+ A^0+ 

since H(r~(si), s^) < H(T+(si), s^) = in view of {in) and (i) in Proposition l5.5l 
In the same way, we get 

lim T(sl + Aa) = lim T{sl + Aa) 

A-!-Ao(a) A^Ao(a) 

= t"(sl + Ao(a)a) < lim T(sL+Aa), 

A^Ao(a) 
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since 'H(r^(s),s) = < 7^(t+(s),s) for all s e C, in view of Proposition 15.51 To 
conclude, we have to establish the smoothness properties put forward in Proposition 
15.61 In view of the monotonicity properties of r i— )■ 'H(t, s) we have just established 
for aU s G /C, aU the three maps are obviously in C^(/C,E;^) nC^{lC,R*^) thanks to 
the implicit function theorem. This concludes the proof of Proposition [5?6l D 
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